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Abstract. We introduce and study the reverse voter model, a dynamics for spin variables similar to 
the well-known voter dynamics. The difference is in the way neighbors influence each other: once 
a node is selected and one among its neighbors chosen, the neighbor is made equal to the selected 
node, while in the usual voter dynamics the update goes in the opposite direction. The reverse voter 
dynamics is studied analytically, showing that on networks with degree distribution decaying as 
k~ v , the time to reach consensus is linear in the system size N for all v > 2. The consensus time 
for link-update voter dynamics is computed as well. We verify the results numerically on a class of 
uncorrelated scale-free graphs. 



INTRODUCTION 

Ordering dynamics is a classical topic in statistical physics, actively investigated for over 
two decades [1]. The best known realization of an ordering process is the dynamics of 
a ferromagnet suddenly quenched from a high-temperature state to a temperature below 
the critical one: ordered domains are formed and coarsen in time until the equilibrium 
state at temperature T is reached. 

More recently, ordering processes have started to be investigated in the completely 
different context of social sciences: the spreading of rumors or beliefs, the dynamics of 
opinions, the diffusion of cultural traits are all phenomena with the common feature 
that, starting from a disordered initial configuration, the dynamics tends to increase 
the similarity of interacting agents, leading toward a more ordered state. Simple (often 
oversimplified) models for such phenomena are akin to models studied in traditional 
statistical physics, but from the point of view of social sciences it is more interesting to 
study their behavior on networks, rather than on regular lattices. 

The voter model is one of the simplest ordering processes [2]. Each individual is fully 
specified by a single 'spin' variable that can assume only two states, say —1 and 1. At 
each time step an individual is chosen at random and his variable is set equal to the 
variable of a randomly chosen nearest neighbor. This step is then iterated many times. 
There are two absorbing configurations for the dynamics: all variables equal to —1 or 
all equal to 1 . Once these consensus states are reached no further evolution is possible. 
When N, the number of agents in the system, is finite, the dynamics reaches sooner or 
later one of the two absorbing configurations. An interesting quantity is then t(N), the 
mean time needed to arrive at full consensus. For regular lattices in d dimensions, this 
time scales as N 2 in d = 1, Nln(N) in d = 2 and as N for d > 2 lHH. 



Some activity has been recently devoted to the study of the voter model on graphs. 
In particular, it has been shown that the time to reach consensus scales as N on the 
Watts-Strogatz graph ]4|,|5|]. On the Barabasi- Albert (BA) scale-free network it is found 
numerically that T grows with /V with an exponent close to 0.88 1 6, 7]. At odds with what 
occurs on regular lattices, voter dynamics does not conserve the average magnetization 
when the number of nearest neighbors (degree) changes from node to node JalS]. In a 
very recent contribution ||9|], Sood and Redner have evaluated analytically the consensus 
time for a generic uncorrelated heterogeneous graph with a scale-free degree distribution 
rik ~ k~ v . It turns out that x(N) scales as N for v > 3, it is a constant for v < 2, and it 
scales as Af(2v-4)/(v-i) f or 2 < v < 3. The presence of a logarithmic correction for v = 3, 
t(/V) ~ N / \n(N) accounts well for the observed exponent 0.88 on the BA network. 

In the elementary step of the original voter model, the selected node A picks at random 
a nearest neighbor B and becomes equal to it. When one considers simple modifications 
of this dynamics there are two possibilities that immediately come to mind. The first, 
similar to the so called Ochrombel simplification of the Sznajd model fllOfl . consists in 
reversing the direction of the update, so that it is the neighbor B that becomes equal 
to the selected node A. We call this dynamics 'reverse' voter model. The second is the 
'link-update' voter model [6]: instead of a site, a link is randomly chosen and the update 
occurs in random direction along it. The link-update dynamics conserves the average 
magnetization for any degree distribution [ 6]. It is clear instead that the 'reverse' voter 
model violates this conservation. 

Whenever the degree distribution is a delta function the three models perfectly coin- 
cide, but this is no longer true on a generic graph. This is evident for example on a star 
graph, where all links connect a hub with the other N — 1 nodes. In a single elementary 
step, the probability that the hub is updated is l/N for the usual voter model, 1/2 for the 
link-update dynamics and 1 — l/N for the reverse voter model. 

In this paper we investigate for which topology the three types of voter dynamics give 
similar results and when, instead, the ordering processes are qualitatively different. 



We consider the reverse voter dynamics on a generic heterogeneous uncorrelated graph 
with normalized degree distribution n\. Let us define p^, the fraction of up spins on 

nodes of degree k, and compute the probability P(k; ► +) that a spin down on a node 

of degree k flips in an elementary update. 

For this flip to occur it is necessary that the spin on the selected node (A) is up (and this 
happens with probability p = Y,k n kPk) an d that the neighbor chosen (B) has degree k and 
is down. The degree distribution of the neighbors of a random node in an uncorrelated 
network is kn k j\i\, where /i ; = Y.k^ n k- The probability that a node of degree k is down 
is 1 — p k , hence 



ANALYTICAL RESULTS 
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Using these expressions one can apply the analytical treatment of Ref. to compute 
the consensus time x(N). The idea is to write down a recursion formula for x as a 
function of the p k . Expanding to second order one obtains 

L^(p,-p)4T-^I^(p + p,-2pp,)^T=l, (2) 
where d k = d/dp k . 

To analyze the role of the different terms in Eq. © we have to consider the dynamics 
of the densities p k . In an elementary update, the number N? p of up spins of degree k, 
changes by an amount dN^ p = P(k; ► +) —P(k; H — ► — ), so that 

dN up 1 k 

* = a& = 7^1,™ - ^ +) - P{k + — >1 = £<P "ft>- < 3 > 

The equation of motion for the magnetization p is 

P=P \n k kp k . (4) 

Mi k 

Starting from their initial values, the densities p k rapidly converge to a stationary value 

p k = p, Vfc. (5) 

This stationary value can be related to the quantity R = /ii E/t tP^ tnat ls conserved 
by the dynamics 

k k c 

where C = 

It is now possible to rewrite Eq. © in the stationary state. Using Eqs. © and © the 
first term is seen to vanish; one is left with 

R(l-CR)dlx = -N. (7) 

The integration of Eq. Q, with boundary conditions t(7? = 0) = and x(R = 1/C) = 0, 
yields 



x(N) = -N 



(8) 



R]n{CR) + -RJ ln(l -CR) 
For an uncorrelated initial configuration, i.e. p k (0) = p(0) 

T(7Y) = -iV/ii/i-i [p(0)lnp(0) + (l -p(0))ln(l -p(0))]. (9) 

Since the prefactor is finite the time to consensus diverges linearly with N — > 00 
for all v > 2. For v = 2 a logarithmic correction is present: x(N) N\n(N), while 
T(iV) -/V 1 /^- 1 ) for 1< v < 2. 
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FIGURE 1. Main: Double logarithmic plot of the consensus time z(N) as a function of N for the reverse 
voter dynamics on a growing network with redirection. The thick line has slope 1 . Inset: The same data 
divided by Nflifl-i. The thick line is here the result predicted by Eq. (|9}. 



Applying the same analytical approach, the consensus time for link-update dynam- 
ics Ja] can be computed as well. The probability of flipping a down spin on a node of 

degree k has the form P(k; ► +) = (Q-tA(1 — pk), where CO = Y,k n kkpk- The consensus 

time turns out to be 

T(iV) = -iV[(l-p)ln(l-p)+plnp]. (10) 

Also in this case the time to consensus depends linearly on the system size N. Remark- 
ably, there is no dependence at all on the network parameters. 



NUMERICAL RESULTS 

In order to validate the anaytical results we have performed numerical simulations of the 
reverse voter dynamics on two classes of scale-free graphs. 

In Fig. [T] we present the consensus time for reverse voter dynamics on a growing 



network with redirection [11] for several values of the exponent v of the scale-free 
degree distribution. 

From the main part of the figure it is clear that, while the asymptotic behavior is, as 
expected, proportional to the system size, sizeable corrections occur for smaller values 
of N. This might be caused by the residual preasymptotic dependence of /ii and /i i on 
N, for the relatively small values of N considered. In the inset we have therefore plotted 
t(N)/(NhiH-i) and compared to the value ln(2) = 0.6931..., predicted by Eq. ©. The 
agreement is not satisfactory. This suggests that correlations present in the network have 
an important influence on the ordering process; they do not alter the asymptotic linear 
dependence on N but affect rather strongly the prefactor. 
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FIGURE 2. Main: Double logarithmic plot of the consensus time t(N) as a function of N for the reverse 
voter dynamics on the uncorrected scale-free networks of Ref. 1 12],with minimum degree equal to 3, for 
Pifc(O) = P(0) = 1 /2- The thick line has slope 1. Inset: The same data divided by jUijU i . Here the thick 
line is the result predicted by Eq. 0, ln(2). 



We have then considered the evolution of the reverse voter model on another class of 
graphs, recently introduced by Catanzaro et al. 11211 . Such graphs are built according to 
the usual algorithm for the configuration model, with the additional prescription that no 
node can have a degree larger than N l l 2 . This modification has the goal of avoiding the 
correlations of the node connectivities present in the usual configuration model. 

It is evident from the results, displayed in Fig. El that even on this class of networks 
one does not see a clean linear dependence on N for the system sizes considered here. 
The effective exponent for small N is larger than one; a crossover to the expected asymp- 
totic behavior is seen for larger values of N. The plot of t(N) / (Af/li/l-i), presented in 
the inset, indicates that the bending in the main plot actually comes from the residual 
Af-dependence of jl\ and jit_i. The quantity t(N) /(Njli/J—i) goes asimptotically to a 
constant value that, as predicted from Eq. © does not depend on the exponent V. The 
only feature of Fig. El not in agreement with the theory is the value of the constant, that 
is actually of the order of two times the expected value, ln(2). 

On the uncorrelated scale-free graphs we have also studied numerically the ordering 
of the link-update voter model. Although Eq. (TTOb predicts the consensus time to be 
completely independent from v we find that the curves do not superimpose (Fig. |3l). 
A simple linear fit yields exponents ranging from 0.96 to 1.02, but a more careful 
inspection indicates that the curves are slightly bent. 

Further insight is provided by the inset, where the consensus time divided by the 
system size N is plotted as a function of N. A transient regime, for small N, with non 
constant t(N) /N is followed by plateau, indicating that the asymptotic linear regime is 
reached. At odds with Eq. (flTJl) . the height of the plateau depends on the value of v, 
ranging from approximately one to two times the expected value ln(2). 

Hence we find that also in this case the analytical agreement correctly captures the 
linear behavior independent from v, while it fails for what concerns the prefactor. We 
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FIGURE 3. Main: Double logarithmic plot of the consensus time t(AQ as a function of N for the link- 
update voter dynamics on the uncorrelated scale-free networks of Ref. II 211 with minimum degree 3, for 
Pfc(O) = P (0) = 1/2. Inset: The same data divided by N. The thick line is the result predicted by Eq. (1101 . 
ln(2). 



do not have a clear understanding of the origin of this disagreement. One possible 
reason could be the presence of residual correlations in the modified configuration model 
that we have used. However, the mismatch occurs also for values of v so large that 
correlations are surely negligible. It is more likely that the different prefactors are an 
effect of the term containing the first derivative in Eq. 0, that we have neglected in the 
analytical treatment, but that may play a role in the initial stage of the ordering process. 



CONCLUSIONS 

We have studied analytically and numerically the ordering behavior of the reverse and 
the link-update voter models, finding that the time to reach consensus grows linearly 
with the system size for scale-free networks with v > 2. 

It is interesting to notice that, for the reverse voter model, the time to reach consensus 
gets larger as v is decreased, since }l\ grows. This can be qualitatively understood as 
follows: in the reverse voter dynamics it is very easy that a hub changes its state, being 
chosen as the neighbor of the selected node. In the usual voter dynamics, instead, the 
presence of a hub favors order, since its many neighbors tend to become equal to it. 
Correspondingly, the time to consensus in the normal voter dynamics is reduced as v is 
decreased 0]. 

As mentioned above, the three types of voter model are perfectly equal when the 
degree is the same for all nodes. This can be seen also from Eqs. © and (ITOb . For 
v — > °°, Hi = (/ii)' and in all cases (including the usual voter dynamics), the consensus 
time is given by Eq. (fill . For finite values of v, instead, the usual voter dynamics leads 
to consensus most quickly, and the reverse dynamics is the slowest (Fig.|4|). 

All the results presented here were obtained under the assumption that the dynamics 
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FIGURE 4. Main: Plot of the consensus time T as a function of v for the three types of voter dynamics 
on the uncorrelated scale-free networks of Ref. fl^l with minimum degree 3, for Pk(0) = p (0) = 1 /2. The 
upper curves are for N = 3000, the lower curves are for N = 300. 



takes place on an uncorrelated network. An interesting issue for future work is the 
investigation of voter dynamics on correlated graphs. 
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